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We develop a quantum process tomography method, which variationally reconstruct the map of a 
process, using noisy and incomplete information about the dynamics. The new method encompasses 
the most common quantum process tomography schemes. It is based on the variational quantum 
tomography method (VQT) proposed by Maciel et al. in arXiv:1001.1793[quant-ph]. 
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The characterization of a quantum system and its dy- 
namics is a daunting challenge. The first question which 
arises in this scenario is what information should we pos- 
sess to characterize the dynamics. To answer this ques- 
tion, one needs to choose a quantum process tomography 
(QPT) scheme. Each procedure demands different re- 
sources and operations. 

There are four general types of QPT procedures: (i) 
standard quantum process tomography (SQPT)[1]; (ii) 
ancilla- assisted process tomography (AAPT)[2-4]; (iii) 
direct characterization of quantum dynamics (DCQD)[5- 
7]; (iv) selective and efficient quantum process tomogra- 
phy {SEQPT)[8-10]. 

In (i) the information is obtained indirectly, perform- 
ing a set of quantum state tomographies {QST)[11— 16] of 
the linear independent states, which spans the Hilbert- 
Schmidt space of interest, after the action of the unknown 
map. The second scheme (ii) - also an indirect procedure 
- makes use of an auxiliary system. The information is 
then extracted by means of QST of the joint space (sys- 
tem and ancilla) . The third one (iii) obtains the dynami- 
cal information directly - by means of quantum error de- 
tection (QED)[17] concepts - measuring stabilizers and 
normalizers. Finally, the last method (iv) - which also 
measures the parameters directly - consists in estimating 
averages over the entire Hilbert space of products of ex- 
pectation values of two operators. For the special case of 
one-parameter quantum channels, there is also an inter- 
esting method developed by Sarovar and Milburn [18]. 

In [15], we developed a variational quantum tomogra- 
phy method (VQT), which has the virtue of yielding a 
state which is a universal lower bound estimator for every 
observable. The VQT approach can reconstruct the state 
with high fidelity out of incomplete and noisy informa- 
tion. The method was successfully employed in a quan- 
tum optics experiment, where entangled qutrits were gen- 
erated [16]. In this letter, we extend VQT to the tomog- 
raphy of quantum processes. The new method inherits 
all the advantages of the VQT, and opens the door to the 
characterization of maps in larger systems, where myr- 
iad of measurements could be necessary. With the varia- 



tional quantum process tomography method (VQPT) we 
propose here, maps can be reconstructed, with high fi- 
delity, using just a fraction of the effort employed in an 
Informationally Complete Measurement. 

The method we derive has the particular form of a lin- 
ear convex optimization problem, known as Semidefinite 
Program (SDP), for which efficient and stable algorithms 
are available [19-21]. SDP consists of minimizing a lin- 
ear objective under a linear matrix inequality constraint, 
precisely. 



minimize c' x 



subject to < F{x) = Fq 



Y^ x,F, > 0, 



(1) 



where c e C™ and the Herniitian matrices Fi e (jnxn 
are given, and x e C™ is the vector of optimization vari- 
ables. F{x) > means that F{x) is a positive matrix. 
The problem defined in Eq.l has no local minima. When 
the unique minimum of this problem cannot be found an- 
alytically, one can resort to powerful algorithms that re- 
turn the exact answer [20]. To solve the problem in Eq.l 
could be compared to finding the eigenvalues of a Her- 
mitian matrix. If the matrix is small enough or has very 
high symmetry, one can easily determine its eigenvalues 
on the back of an envelope, but in other cases some nu- 
merical algorithm is needed. Anyway, one never doubts 
that the eigenvalues of such a matrix can be determined 
exactly. 

A bona fide completely positive and non increasing 
trace map £ can be generally represented as 



£ip) 



E 



X^lE,pE], 



(2) 



where p is the system initial state and the {Em} form an 
IC-POVM (Informationally Complete Positive Operator 
Valued Measure), i.e. a complete basis in the Hilbert- 
Schmidt space satisfying 



i=l 



(3) 



The {xij} defines the super-operator x, which has all the 
information about the process. It is a Hermitian positive 
operator. Thus the super-operator can be thought as a 
cP X cP density matrix in the Hilbert-Schmidt space with 
df^ independent real parameters (or d"^ — d^ in the trace 
preserving case). 

Now we recast both SQPT(i) and AAPT(h) - which 
rely on tomography of states - using the VQT[15] 
methodology. We will name the output states of the un- 
known map as g'^ = £{p^) 



Ylti=iXijErP^E]. Suppose 



elements of the POVM (n'"' < d^) in the fc*'' output 



state have been measured, namely 



TriExp'^) 



pI 



Xe [l.n'']. 



(4) 



Note that p^ are positive numbers, the known probabili- 
ties. The probabilities obtained from an experiment are 
noisy, thus: 

(1 - A1)p1 < TritEl) < (1 + A\)pI a G [l,n\ (5) 

with A| positive and hopefully small. Let us refer to the 
space spanned by the unmeasured POVM elements as the 
unknown subspace. Then we can define the Hamiltonian 



^'= E ^. 



(6) 



Thus we obtain a cost function over the unknown sub- 

space, namely, Tr{g''H'') = Tr{ ^ XijE,p^E]H''). 

This linear functional should be minimized, for we do not 
know the action of the map on the unknown subspace. 

The SQPT method demands quantum state tomog- 
raphy in all linearly independent states which span the 
Hilbert-Schmidt space. With n measurements in kt dif- 
ferent states, the variational SQPT reads: 



minimize 



{Y,Tr{~Q^H'^) 






x>o, 

Tr(e^) < 1, 



subject to < 



A*;>o, 



VA'^ G [1 



A^')p^ < Tr{g''E'l) < (1 + A^)p 



\)f\-' 



and /c — [1, kt] 



(7) 



Eq.7 returns a map x which is the optimal approximation 
to the unknown process x- Note that, at the same time, 
we were able to identify g^ optimally. In Fig.l, we illus- 
trate the aplication of the method for the reconstruction 
of two-qubit processes. 

In the AAPT method, one adds an ancillary system 
with the same dimension of the main one. Then the 
quantum process takes place in half subspace, and finally 
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FIG. 1: Reconstruction of two-qubit processess of all ranks. 
We plot the number of independent POVM elements neces- 
sary to reconstruct the map against the process rank. For 
each rank, we randomly generated 100 processes. 



a quantum state tomography is performed in the whole 
space, ancilla plus main system. The output state now 



reads g = {K^ £){p) = Y^ij^^ X»j(I ® E,)p^{V ® Ej)l 
With n measurements performed in this scheme, AAPT 
can be recast as 

n 

minimize {Tr{gH) + \^ A^) 

A=l 

X>0, 
Trig) < 1, 
subject to ^ Ax > 0, 

(1 - Aa)pa < Tr{gEx) < (1 + Aa)pa, 
VAe [l,n]. 

(8) 

In conclusion, we have introduced a new method 
to perform quantum process tomography, which recon- 
structs a map, with high fidelity, using noisy and incom- 
plete information. The method is linear and convex, and 
its unique solution can be obtained very efficiently. It 
opens the door to the characterization of the dynamics of 
larger quantum systems, avoiding the need of very large 
informationally complete measurements. 
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